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Abstract 

The Higgs potential of the standard model with an additional real Higgs singlet is studied 
in order to examine if it may allow the strongly first order electroweak phase transition. It 
is found that there are parameter values for which this model at the one-loop level with 
finite temperature effect may allow the desired phase transition. Those parameter values 
also predict that the masses of the neutral scalar Higgs bosons of the model are consistent 
with the present experimental bound, and that their productions in e + e~ collisions may be 
searched at the proposed ILC with yfs — 500 GeV in the near future. 



I. Introduction 



The possibility of baryogenesis by means of electroweak phase transition has recently been 
widely examined, since the electroweak baryogenesis can in principle be tested in the future 
accelerator experiments [1]. If the electroweak phase transition is strongly first order, it can 
fulfil the departure from thermal equilibrium which is one of the three conditions required by 
Sakharov that are necessary for the dynamic generation of the baryon asymmetry during the 
evolution of the universe [2]. It has already been observed that in the standard model (SM) the 
electroweak phase transition cannot be strongly first order unless the mass of the scalar Higgs 
boson is smaller than its lower bound set experimentally by LEP [3]. The sufficient strength 
of the first order electroweak phase transition is essential for preserving the generated baryon 
asymmetry at the electroweak scale. In the literature, a number of articles have been devoted 
to study the possibility of accommodating the strongly first order electroweak phase transition 
in various models beyond the SM [4]. 

Among them, an interesting possibility has been investigated several years ago, where an 
extension of the SM with a real Higgs singlet field has been adopted within the context of the 
electroweak phase transition [5]. We consider that model is inspiring, because adding a real 
Higgs singlet field is the simplest extension of the Higgs sector of the SM. In that model, the 
strength of the first order electroweak phase transition has been stronger than that in the case 
of the SM. Due to the presence of cubic terms in the tree-level Higgs potential, a strongly first 
order electroweak phase transition has been ensured for a scalar Higgs boson mass of around 
100 GeV, for the top quark mass of 127 GeV. 

In this paper, we would like to reexamine the model in a more rigorous way, in the sense 
that we perform the full integration of the finite-temperature effective potential instead of ap- 
proximating it up to second order in temperature, that we include the Higgs contributions to 
the finite temperature Higgs potential at the one-loop level, and that we take into account the 
zero-temperature one-loop Higgs potential. As the experimental lower bound on the scalar Higgs 
boson mass by now is improved to be more than 114 GeV, and the top quark mass is usually set 
as about 175 GeV, we would like to explore the parameter space of the model, with these new 
numbers in mind, where the first order electroweak phase transition may be strong enough. We 
find that the model is able to produce strongly first order electroweak phase transition for some 
parameter regions where a scalar Higgs boson may have a mass larger than 114 GeV. Further, 
we calculate cross sections for scalar Higgs productions in e + e~ collisions. 

II. The Higgs Potential 

Let us consider the Higgs potential proposed by Ref. [5] where one real Higgs singlet field 
is introduced to the SM Higgs sector. We follow the notations of Ref. [5]. Thus, the Higgs 
potential consists of one complex Higgs doublet field 4> and one real Higgs singlet field S. The 
tree-level Higgs potential is given by 

V°(4>, S) = X^U) 2 - fifa + ^ - f - ^S 3 + 2A(^)S 2 - f , (1) 

where a, co, /j,^, and us are parameters with mass dimension, and the rest are dimensionless. 
We assume that A^A^, and A have a positive values. The real part of the neutral component 
of the Higgs doublet field is, as in the SM, going to develop the non-zero vacuum expectation 
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value (VEV) u = V2(4>) while the Higgs singlet field v = (S). The tree-level values for u and v 
are determined by the extremum conditions 

dV°{u,v) _ dV°{u,v) _ 



du 



dv 



where 
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We denote the zero temperature values of u and v as uq and vq, respectively. We take uq = 246 
GeV whereas vq is not fixed but variable. 

In this model, there are two neutral scalar Higgs bosons. At the tree level, the masses for 
them are given by 
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We associate the lighter neutral scalar Higgs boson with the real part of the neutral component 
of the Higgs doublet <j) while the heavier one with the real Higgs singlet S. 

At the tree level, the gauge bosons and the top quark mass are mw{uo) = 80.423 GeV, 
m-ziuo) = 91.187 GeV, and m t (uo) = 175 GeV. They are all defined in terms of uo, independently 
of -uo- 

Now, we consider the one-loop contributions. At the one loop level, the zero temperature 
effective potential [6] and the finite temperature one [7] are respectively given by 



64tt 2 
2vr 2 



log- 



mf(u, v) 



A 2 



dx x 2 log 



1 ±exp (~\Jx 2 + mf(u,v)/T 2 



(6) 
(7) 



where nw = 6, nz = 3, rit = —12, = ns = 1 account for the degrees of freedom of each 
participating particles. The negative sign in V T is for bosons while the positive sign for fermions. 

In V 1 , we take the renormalization scale as A = 246 GeV. The A dependence of the potential 
would disappear if all of higher-order contributions be taken into account. This is desirable 
since physical observable quantities should be independent of the renormalization scale. As we 
consider only the one-loop contribution in V 1 , the scale dependence remains. 

In order to reduce the effect of the scale dependence, one may introduce appropriate counter 
terms to the one-loop contribution, after regularizing the relevant fields in the effective poten- 
tial. In our model, the scale dependence on the mass of top quark or on the mass of weak gauge 
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bosons can completely be eliminated in the one-loop contribution. This is possible because 
the Higgs singlet in our model does not contribute to fermion masses or the weak gauge boson 
masses. Thus, the one-loop contributions to these masses in our model are exactly identical to 
those in the SM. Hence, the scale dependence of the one-loop contributions to these masses may 
disappear if mf(u)log(A 2 ) in Eq. (6) is substituted by mf(u) log(m|(uo)) — 2m 2 (u)m 2 (uo) « 
mj(u) log(m|(uo)), where we can neglect 2m 2 (u)mf (uo) . This is effectively equivalent to replac- 
ing A by nii(uo). 

The case of the Higgs boson mass is more complicated than the case of the fermion masses 
or the weak boson masses, since the squared mass of the Higgs boson at the tree level depends 
not only u but also v, and it is given by Eq. (5) which is not in the form of C\ + C2U 2 (where 
C\ and C2 are constants). Nevertheless, comparing with the one-loop contribution to the top 
quark mass, we might effectively set A in the one-loop contribution to the Higgs boson mass by 
either m^(uo, vq) or m° s (uQ, vq), in a similar way as A — > m t (uo) in the case of top quark mass. 

If a model is related to a new physics at certain scale, the renormalization scale might be 
taken comparable to the energy scale of the new physics, in order to consider the higher-order 
effects of the model reasonably. For example, in a supersymmetric axion model, A is taken to be 
10 11 GeV in the one-loop effective potential, where new physics is expected to emerge at that 
energy scale [8]. In our study, we are interested in the tree-level mass of the Higgs boson at the 
electroweak symmetry breaking scale. Thus, in our case, we might set A = 246 GeV. 

For small rrii/T, that is, at high temperature, V T can be expanded into the well-known form 
for the temperature induced effective potential. Rather, we perform the integration without 
using the high temperature approximation in our analysis, in order to consider relatively large 
scalar Higgs boson mass. 

The full effective potential that we consider therefore consists of 

V(u, v) = V°(u, v) + V 1 ^, v) + V T (u, v) . (8) 

Note that the effective Higgs potential possesses a symmetry under interchange u <-> —u. Thus, 
we may confine ourselves within the region of u > 0. 

In order for the first order electroweak phase transition to take place, we need two separate 
vacua. The vacua are defined as the points in the (u, v)-plane at which the Higgs potential has 
local minima. The two extremum conditions for the Higgs potential may be expressed as 

dV(u,v) dV(u,v) 

— — — = fo{u,v) = ufi(u,v) = , — — — = f 2 ( u ,v)=0. (9) 

where the first equation is decomposed such that it is explicitly satisfied by u = 0. 
In the above expressions, f±(u,v) and f 2 {u,v) are given as 
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/(a, 6) = {a log A _ ft log ^} + 1 , (11) 
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(12) 



Note that the negative sign in df(m 2 ) is for bosons whereas the positive sign is for fermions. 

Since the first extremum equation is satisfied by either u = or fi(u,v) = 0, one of the 
two points in the (u, v)-plane which satisfy the two extremum equations is simply (0, vi) where 
v\ is the solution of /2(C), v\) = 0. This is the vacuum of the unbroken phase state. We solve 
72(0,^1) = to obtain v\ numerically by the bisection method. Let the other point in the (u, 
-u)-plane which satisfies the two extremum equations be denoted as (u 2 , v 2 ), where u 2 7^ 0. This 
is the vacuum of the broken phase state. This point is obtained by the Newtons's method. 

In order to check that the extremum (u 2 ,v 2 ) thus obtained is minimum indeed, we need a 
Jacobian matrix 

dfi{u,v) dfi{u,v) 

(13) 



Jij(u,v) = 



du ' dv 
df 2 (u,v) df 2 (u,v) 



du dv 

This point is a minimum of the Higgs potential if they satisfy Ju(i = 1, 2) > and det{Ji 3j 



> 0. 



The distance between the two minima in the (u, v)-plane is defined as v c = \f{u2) 2 + (v 2 — ^l) 2 - 
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The nonzero v c ensures that the electroweak phase transition is discontinuous, that is, first 
order. The strength of the first order electroweak phase transition is measured by the critical 
temperature T c , which is defined as the temperature at which the values of the potential at the 
two minima are equal: V(0,v±) = V(u2,vz). We determine T c by varying temperature T until 
V(0, v\) = V(ii2, V2) is met. If the ratio v c /T c is larger than 1, the first order electroweak phase 
transition is defined conventionally as a strong one. 

Now, let us evaluate the masses of the two neutral scalar Higgs bosons from the symmetric 
2x2 mass matrix at the one-loop level. The masses are evaluated at zero temperature. Thus, 
V T does not contribute to the masses at the one-loop level, and we should take u = uq and 
v = v in the final expressions. The elements of the mass matrix for the neutral scalar Higgs 
bosons are given by 



M22 = (M^ + M^) 
M 12 = (M° 12 + Ml 2 ) 



(14) 



where and are obtained from V° and V 1 , respectively. Explicitly, they are obtained as 
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where 



g(a,b) = ^±\log- + 2, (17) 
a — b a 



A 3 = 2 | (3A S - \)v + j - a j + ^2A^u 2 - 2X s v 2 - + aw 




u 2 



A 4 = 



- 2(3A - A) |(3A 5 - X)v + j - a J u + 8A (^2At; - ^ u . (18) 



In terms of these matrix elements, the neutral scalar Higgs boson masses at the one-loop 
level are given as 



m % m s = \ 



(Afn + M 22 ) T ^/(M 22 - Mn) 2 + (M 12 ) 2 



(19) 



which depend on the quartic coupling coefficients A^, Xs, and A, and three parameters vq, w, and 
a, having mass dimension. Here, it is assumed that ms < ms- The quartic coupling coefficients 
satisfy the condition of Xi/4ir < 1 (Aj = X^, Xs, X) for perturbation theory to work, in principle. 
Thus, they are not large. In practice, we may set their values to be order of 1 or less. 

The value of vo should be comparable to that of uq = 246 GeV, because it has a similar 
property as uq. Note that u is the Higgs field in the SM for v = and uo in the SM is the 
electroweak symmetry breaking scale. Likewise, we assume that the other parameters u> and a 
are not so much different from uq. Since we are interested in the mass region of 200 to 400 GeV 
for the upper bound on the Higgs boson masses, we choose the values of the parameters such 
that they yield the Higgs boson masses within this region. 

For numerical study, we randomly search the parameter space by varying the relevant pa- 
rameters within the ranges of < A^,As,A < 0.7, < vq < 300 GeV, and < — LO,a < 100 
GeV. Starting with a set of relevant parameters whose values are chosen randomly within the 
allowed ranges, we examine if v c /T c > 1, which is the criterion for the first order electroweak 
phase transition to be strong. If this criterion is met, we continue to calculate the masses of 
the neutral scalar Higgs bosons. This sequence of calculations produces a point in the (m^, 
v c /T c )-plane or in the (ms, f c /T c )-plane. In this way, we perform calculations to obtain 5000 
points. Figs, la and lb show the results of our numerical studies on the masses of the neutral 
scalar Higgs bosons. In Fig. la, the 5000 points thus obtained are plotted in the (ms, v c /T c )- 
plane, and in Fig. lb in the (ms, t> c /T c )-plane. Notice that the points in the (ms, i> c /T c )-plane 
are scattered slightly toward to the right of the plane, as compared to the points scattered in 
the (m^, -u c /T c )-plane, since we assume that ms < ms- These figures, therefore, indicate that 
there are sufficiently wide regions in the parameter space where the first order electroweak phase 
transition is strong as well as the masses of the neutral scalar Higgs bosons are relatively large. 

We now examine the possibility of detecting these neutral scalar Higgs bosons in the future 
linear e + e~ collider. The three main production processes for the neutral scalar Higgs bosons 
in e + e _ collisions are: 

Higgs-strahlung : e + e~ — ► Zcf>, ZS , 
WW fusion : e + e~ — ► v e v e (\>, v e v e S , 
ZZ fusion : e + e~ — > e + e~d>, e + e~ S , 
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where the Higgs-strahlung process is dominant at the center of mass energy corresponding to 
that of LEP2 whereas WW and ZZ fusion processes are comparable at much higher center of 
mass energy, such as that of the future linear e + e~ collider with > 500 GeV. 

The cross section for 4> or S production via the Higgs-strahlung process is related to the 
corresponding production cross section of the SM scalar Higgs boson, H, as 

cr(Z(j)) = cos 2 (-f)a S M(ZH) , 

a(ZS) =sm 2 ( 1 )a SM (ZH) , (20) 

and, similarly, the cross sections for 4> and S production via the WW fusion process or the ZZ 
fusion process are related to the corresponding production cross section of the SM scalar Higgs 
boson, H, as 

a(VV<P) = cos 2 7 a SM {VVH) , 

a(VVS) = sin 2 7 g S m{VVH) , (21) 

where VV = WW or ZZ, and 7 is the mixing angle between the two neutral scalar Higgs 
bosons: 

o M n - M 22 

cos 27 = . (22) 

^{M 22 - M xl y + A{M l2 y 

These relationships hold because the coefficients for VV<j> and VVS coupling in our model are 
respectively cos 7 and sin 7 times the VVH coupling coefficient in the SM. 

For a given set of parameter values in our model, the above relationships indicate that, in 
case of the Higgs-strahlung process, although neither o{Z<$) nor a(ZS) is larger than cfsm(ZH), 
they cannot be simultaneously smaller than the SM value. If one becomes negligible, the other 
approaches to the SM value. The same argument holds in case of the fusion processes. If 
cos 7 = 1, the neutral scalar boson <j> in our model would behave like H in the SM. Its properties 
would be identical to those for H in the SM. In particular, if cos 7 = 1, the lower bound on the 
mass of <j> would be the same as the lower bound on the SM Higgs boson mass, 115 GeV, which 
has been set by the LEP experiments. If, on the other hand, cos 7 = 0, the neutral Higgs boson 
S in our model would behave like H in the SM. 

Let us consider the Higgs-strahlung process. At the center of mass energy of LEP2, the 
Higgs-strahlung process in the SM is dominant over the fusion processes. We first calculate 
both cr(Z(f>) and a(ZS), and then choose the larger one 

o"max = max{a(Z<f>),(r(ZS)} , 

since at least one of them is comparable to <tsm(ZH). We calculate cr^ ax for the each of the 
5000 sets of relevant parameter values which produce the points in Figs. 1. We plot the result 
against in Fig. 2. The number of points in Fig. 2 is smaller than those in Figs. 1, because 
those points in Figs. 1 yielding > 118 GeV are kinematically excluded from Fig. 2. 

From Fig. 2, the lower bound on might be roughly estimated. It can be seen from Fig. 
2 that the lowest value for m<f, is about 20 GeV. However, this lower bound on m<f, does not 
take into account the discovery limit. If the discovery limit of the neutral scalar Higgs boson 
4> in our model via the Higgs-strahlung process in e + e~ collisions at the LEP2 energy is, for 
example, 400 fb, those points above 400 fb should be excluded since LEP2 had not discovered 
(j) or S. Then, the lower bound on is pushed upward to about 65 GeV. It is known that the 
cross section for the SM Higgs production via the Higgs-strahlung process in e + e~ collisions at 
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the LEP2 energy is about 200 to 250 fb for mn ~ 110 GeV. Assuming that the discovery limit 
of the <fr production cross section in our model is roughly the same order of magnitude, namely, 
a few hundred fb, we might estimate that > 60 GeV. Similar argument would apply to the 
lower bound on the mass of the neutral scalar boson S, if we plot <r^ ax against mj. 

The LEP2 experiments have provided detailed, mass dependent, limits on the Higgs produc- 
tion and decay cross sections, yielding a limit on the Higgs mass as a function of the ratio of 
the cross section times the branching ratio of its decay into bb pairs to the SM value. Combined 
with similar analyses for other decay modes, the LEP2 experiments may exclude a Higgs boson 
with mass of order 60 GeV even if its production cross section is a tenth of the SM value. In 
Ref. [9], the same analysis has been done independently of the specific hadronic decay mode 
of the Higgs boson. Comparing the result of Ref. [9] with our values, one can conclude that 
the bound on is indeed of order 60 GeV. Among the points in Fig. 2, those points with 
Tn4> > 60 GeV and small cross section, lower than 100 fb, are consistent with LEP data. How- 
ever, the key point here is that the LEP2 experiments cannot exclude the whole points in Fig. 2. 
This implies that there are parameter regions in our model, where both the electroweak phase 
transition is strongly first order and the productions of the neutral scalar Higgs bosons via the 
Higgs-strahlung process in e + e~ collisions are consistent with the LEP2 experiments. 

Now, we increase of the e + e~ collisions up to the proposed ILC center of mass energy of 
500 GeV. At this energy, depending on the scalar Higgs boson mass, the WW and ZZ fusion 
processes are comparable to the Higgs-strahlung process. At this energy, the search for the SM 
Higgs boson via the Higgs-strahlung process in e + e~ collision may extend to the mass region 
up to about 400 GeV. The discovery limit of the cross section for the SM Higgs production via 
the Higgs-strahlung process at the ILC of yfs = 500 GeV is about 3 fb. For our model, taking 
the mixing angle factor cos 2 7 or sin 2 7 into account, the discovery limit would slightly be larger 
than the SM one. It is thus possible to assume that the discovery limit of the cross section for 
the Higgs production via the Higgs-strahlung process in our model is less than 10 fb. 

At yfs = 500 GeV, all of the 5000 sets of parameter values we select in Figs. 1 are kinemat- 
ically allowed to produce or S in e + e _ collisions. We calculate the cross sections for 4> an d 
S production via the Higgs-strahlung process, the WW process, and the ZZ process, namely, 
a(Z</)), a(ZS), a{WW<f>), a(WWS), a\zZ<f>), and a(ZZS). We also choose in each channel the 
larger one between the (j) production cross section and the S production cross section: c^ax = 
max{a(Z</)),a(ZS)}, a^ ax = m&x{a(WW<f>), a(WWS)}, and a^ ax = m&x{a(ZZ(/)),a(ZZS)}. 

Our results are plotted in Figs. 3-5. Figs. 3 show the results for the Higgs-strahlung process, 
Fig. 3a for a(Z(f)), Fig. 3b for a(ZS), and Fig. 3c for cr^ ax . For simplicity, all of these figures 
are plotted against m^. Fig. 3c suggests that, at y/s = 500 GeV, the production cross sections 
via the Higgs-strahlung process are perceivably higher than the discovery limit of 10 fb, for 5000 
sets of parameter values in our model. Therefore, these sets of parameter values in our model 
may be explored at the ILC with ^fs = 500 GeV. 

Figs. 4 show the results for the WW fusion process, and Figs. 5 for the ZZ fusion process. 
The production cross sections via the ZZ fusion process are generally smaller than those via 
the other channels. The production cross sections via the WW fusion process are considerably 
large for smaller scalar mass, but rapidly decrease as the scalar mass increases. These behaviors 
are expectably consistent with the SM predictions. 
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III. Discussions 



Successful baryogenesis requires the presence of CP violating sources beyond the CKM phase 
in the SM. This is because the CKM phase in the SM cannot produce sufficient CP violation 
for the required amount of the electroweak baryogenesis. The Higgs sector of our model cannot 
accommodate CP violation since it consists of the SM Higgs doublet and a real Higgs singlet. 
Thus, we need sources of CP violation other than its Higgs sector. As possible sources of CP 
violation, scenarios for new physics are good candidate. Assuming the existence of new physics 
have usually been induced to solve the hierarchy problem. 

Among various scenarios for new physics, two of them can be incorporated with our model 
[10, 11]. In one scenario [10], higher-dimensional interaction terms associated with the Higgs 
doublet may be introduced at the electroweak energy scale. In this case, the CP violation 
occurs at the electroweak scale, and the effective Yukawa coupling for top quark is given by 
hf f = ht[l + c t e*{4> 2 - v 2 /2)/A 2 c ], where £ is a CP violating phase, q is a real parameter, and 
A c is a cut off energy scale related to the new physics which is much larger than the electroweak 
energy scale. Ref. [10] suggests that the energy scale for the new physics might be the Majorana 
mass scale of the right-handed neutrino in left-right symmetric models. 

Another scenario is more realistic and suitable for our model [11], where energy scale for the 
electroweak phase transition and the energy scale for the CP violation can be made different. 
The Higgs sector of the model studied in Ref. [11] has a conventional complex Higgs doublet 
and an extra complex Higgs singlet. The conventional complex Higgs doublet play the role of 
the SM Higgs doublet, introducing the SM Higgs boson, and the extra complex Higgs singlet 
introduces two additional neutral Higgs bosons. Both of the additional neutral Higgs bosons are 
heavier than the SM Higgs boson; one of them stays at the electroweak scale whereas the other 
has a mass of high energy scale. Thus, among the three neutral Higgs bosons of the model, 
there are two neutral Higgs bosons at the electroweak scale and the third one is decoupled from 
them. 

The model considered in Ref. [11] introduces not only the extra complex Higgs singlet 
but also extra vector-like down quark. The introduction of a vector-like quarks to the SM is 
attractive since they naturally arise in grand unified theories such as Eq, noted Ref. [11]. In the 
fermion sector of the model considered in Ref. [11], the SM fermions stay at the electroweak scale 
whereas the vector-like down quark stays at high energy scale. Consequently, the electroweak 
phase transition occurs at the electroweak scale, whereas the CP violation occurs at high energy 
scale. In this case, therefore, the contributions of the heavy third neutral Higgs boson and the 
vector-like down quark in the finite-temperature effective potential are Boltzmann suppressed 
and thus negligible when the electroweak phase transition are considered at the electroweak scale. 
In other words, the model considered in Ref. [11] becomes effectively identical to our model at the 
electroweak scale. Actually, by suitably redefining the tree-level Higgs potential, and decoupling 
the third neutral Higgs boson, in the model considered in Ref. [11], it becomes exactly the 
tree-level Higgs potential of our model. Thus, our model for the electroweak phase transition 
may be extended naturally into the model considered in Ref. [10] in order to accommodate the 
CP violation at high energy scale, without modifying our study at the electroweak scale. 

The terms proportional to a and to in the tree-level Higgs potential come from Z2 symmetry 
breaking terms. If a = u = 0, these two terms are absent, and the parameter region allowing 
strongly first order electroweak phase transition disappears. This can be seen, for example, in 
Fig. 2 of Ref. [11], where a is the same and uj is — 4£. Thus, the Z2 breaking terms with 
a 7^ to 7= are crucial for the strongly first order electroweak phase transition. 
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Actually, there should be a potential barrier between two vacua in the tree-level Higgs 
potential in order for the strongly first order electroweak phase transition to occur. As an 
example, in a nonminimal super symmetric model, a trilinear term may be introduced in its tree- 
level Higgs potential to replace the /i-teim. This term allows strongly first order electroweak 
phase transition in the nonminimal supersymmetric model. Also, in Ref. [5], one can see that 
there are a global minimum and a local minimum in the tree-level Higgs potential, which become 
two separate global minima by the contribution of the temperature dependent effective potential. 

IV. Conclusions 

The SM with an additional real Higgs singlet in its Higgs sector is studied within the context of 
the electroweak phase transition. Radiative corrections by the one-loop contribution from gauge 
bosons, top quark, and the Higgs bosons are included, as well as the finite temperature effects. 
The one-loop finite-temperature corrections are integrated numerically, instead of employing 
the high temperature approximation. By randomly choosing a large number of sets of the 
relevant parameter values within the ranges of < A0,As,A < 0.7, < vq < 300 GeV, and 
< -co, a < 100 GeV, we select those sets of parameter values that satisfy v c /T c > 1. For those 
particular sets of parameter values, we calculate the masses of the two neutral scalar Higgs 
bosons of the model and the production cross sections for them via three dominant processes 
in e + e~ collisions. We find that there are as much as 5000 sets of parameter values of this 
model that may accommodate the strongly first order electroweak phase transition, which yield 
reasonably large masses for the neutral scalar Higgs bosons and large production cross sections 
in e + e~ collisions with ^/s = 500 GeV. Thus, for certain parameter values, the prediction on 
the masses of the neutral scalar Higgs bosons are consistent with the present experimental lower 
bound. The results of our calculations suggest that this model may be tested in the near future 
at the proposed ILC. 
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Figure Captions 



Fig. la : Within the ranges of < X^,X S , A < 0.7, < v < 300 GeV, and < -u,a < 100 
GeV, 5000 sets of parameter values are randomly selected which yield v c /T c > 1, and they are 
plotted against m^. 

Fig. lb : For the same 5000 sets of parameter values as in Fig. la, v c /T c against mg are plotted. 

Fig. 2 : For the same 5000 sets of parameter values as in Fig. la, the larger one of the 
production cross sections, cr^ ax = max{cr(Z0), a(ZS)}, for the neutral scalar Higgs bosons via 
Higgs-strahlung process in e + e~ collisions with ^fs = 209 GeV is plotted against m^. 

Fig. 3a : For the same 5000 sets of parameter values as in Fig. la, the cross section for </> 
production via Higgs-strahlung process in e + e~ collisions with ^/s = 500 GeV is plotted against 

771,0. 

Fig. 3b : For the same 5000 sets of parameter values as in Fig. la, the cross section for S 
production via Higgs-strahlung process in e + e~ collisions with y/s = 500 GeV is plotted against 

7770. 

Fig. 3c : For the same 5000 sets of parameter values as in Fig. la, the larger one between the 
cross section for cf> production and the cross section for S production via Higgs-strahlung process 
in e + e~ collisions with yfs = 500 GeV is plotted against m^. 

Figs. 4a-c : The same as Figs. 3a-c, except for WW fusion process instead of the Higgs-strahlung 
process. 

Figs. 5a-c : The same as Figs. 3a-c, except for ZZ fusion process instead of the Higgs-strahlung 
process. 
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Fig. la: Within the ranges of < A^, \ s , A < 0.7, < v < 300 GeV, and < -u,a < 100 
GeV, 5000 sets of parameter values are randomly selected which yield v c /T c > 1, and they are 
plotted against m^. 
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Fig. lb: For the same 5000 sets of parameter values as in Fig. la, v c /T c against ms are plotted. 
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Fig. 2: For the same 5000 sets of parameter values as in Fig. la, the larger one of the 
production cross sections, cr^ ax = max{cr(Z<^), a(ZS)}, for the neutral scalar Higgs bosons via 
Higgs-strahlung process in e + e~ collisions with y/s = 209 GeV is plotted against m^. 
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Fig. 3a: For the same 5000 sets of parameter values as in Fig. la, the cross section for cf> 
production via Higgs-strahlung process in e + e~ collisions with ^/s = 500 GeV is plotted against 

777,0. 
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Fig. 3b: For the same 5000 sets of parameter values as in Fig. la, the cross section for S 
production via Higgs-strahlung process in e + e~ collisions with y/s = 500 GeV is plotted against 

777,0. 
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Fig. 3c: For the same 5000 sets of parameter values as in Fig. la, the larger one between the 
cross section for <fi production and the cross section for S production via Higgs-strahlung process 
in e + e~ collisions with yfs = 500 GeV is plotted against m^. 
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Fig. 4a: The same as Figs. 3a, except for WW fusion process instead of the Higgs-strahlung 
process. 
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Fig. 4b: The same as Figs. 3b except for WW fusion process instead of the Higgs-strahlung 
process. 
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Fig. 4c: The same as Figs. 3c except for WW fusion process instead of the Higgs-strahlung 
process. 
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Fig. 5a: The same as Figs. 3a, except for ZZ fusion process instead of the Higgs-strahlung 
process. 
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Fig. 5b: The same as Figs. 3b, except for ZZ fusion process instead of the Higgs-strahlung 
process. 
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Fig. 5c: The same as Figs. 3c, except for ZZ fusion process instead of the Higgs-strahlung 
process. 
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